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Abstract 

A Born-Markov master equation is used to investigate the decoherence of the state of a macro- 
scopically occupied mode of a cold atom trap due to three-body loss. In the large number limit 
only coherent states remain pure for times longer than the decoherence time: the time it takes for 
just three atoms to be lost from the trap. For large numbers of atoms (N > 10 4 ) the decoherence 
time is found to be much faster than the phase collapse time caused by intra-trap atomic collisions. 
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Decoherence due to environmental coupling JT|, [|, ||] can be very rapid and has the effect of 
making certain pure quantum states unobtainable for all practical purposes. This means that 
the environment largely determines the state of the system as only so called "robust states" 
survive for long periods of time. This is also true for cold trapped atoms and in particular 
when a Bose-Einstein condensate is present in the trap [[|. Even at zero temperature the 
condensate atoms are coupled to the free modes outside the trap by atomic loss processes 
and so these free modes act as an environment for the condensate. It has been argued |5]. |(| 
that environmental coupling will leave the condensate in a coherent state with a well-defined 
phase due to the robustness of the coherent states. 

The form of the environmental coupling is crucial in determining the robust states of an 
open system [|TJ. It is well known that the coherent states are robust under linear coupling 
0,0, or single-body loss, in the sense that an initial coherent state decays in amplitude 
but otherwise remains unchanged. More important for large Bose-Einstein condensates is 



loss due to three-body recombination M, M, 10, 11], which limits the lifetime and size of the 



present condensates. Three-body loss can be considered to be a fundamental loss process 
for trapped atoms as it cannot be reduced (except by changing the scattering length [|12||), 
unlike loss due to spin-relaxation which does not occur in an optical trap and unlike loss due 
to collisions with untrapped atoms which can be reduced by improving the vacuum. In this 
letter I aim to show that three-body loss, gives rise to a very rapid decoherence of certain 
states in spite of the fact that it leads to only a slow decay in the number of trapped atoms. 
In the large number limit the robust states are similar to coherent states with a well-defined 
phase. 

In opposition to the above argument, elastic atom-atom collisions within the condensate 
will give rise to a dephasing between states of different number and in the absence of en- 
vironmental coupling lead to a phase-collapse [JT3[, [14, 15 . However, I will show that the 
timescale at which the three-body loss decoheres an initial state to the robust states is much 
more rapid than the phase collapse for large condensates. The collisions are then more likely 
to lead to a slow diffusion of a well-defined phase rather than a collapse. 

To investigate the quantum properties of three-body loss I derive a Born-Markov master 
equation for the density matrix of the trapped atoms. Born-Markov master equations have 
been an important tool for describing radiative decay in quantum optics (see for example 

g). 
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For cold atoms in a trap the three-body recombination process produces a molecule and 
a single atom, (e.g. Rb + Rb + Rb — > Rb + Rb2), where the molecule is an atom-atom 
bound state. The weakest binding energy, E b , is related to the s-wave scattering length, 
a, by Ef, = hcui, = h 2 /ma 2 ||17|. In this work I consider only the repulsive interaction case, 



a > 0. Due to the large binding energy and momentum conservation the molecule and atom 
are both created with a large kinetic energy and escape from the trap. We are then interested 
in the influence of this loss process on the behavior of the remaining trapped atoms. 

For low-temperature atoms the three-body process becomes a point interaction charac- 
terized by a single parameter k such that the Hamiltonian can be written 



H i = ?ikJ dx [V4(x)^(x)4(x) + h.c] , (1) 



where ^t( x ), ^f( x ) an d 0f( x ) are the field operators for the trapped atoms, the free (un- 
trapped) atoms and free molecules, respectively. The field operators commute with each 
other but have the usual Bose commutation relations with their Hermitian conjugate. 

The Hamiltonian for the cold trapped atoms, if-r, is the usual many-body Hamiltonian 
describing two-body s-wave scattering (see, for example, the review article fl8f ). In most 
circumstances three-body elastic collisions will be negligible in comparison to two-body 
elastic collisions so I do not include them here. 

The density matrix of the total system satisfies the interaction-picture equation of motion 
dptot(t)/dt = -i[iJ / (t),p tot (t)]/^, where H^t) = e -i(HT+H F )t/n Hie i(H T +H v )t/n and Hp ig the 

Hamiltonian of the free particles. The free particles and the trapped atoms are assumed to 
be initially uncoupled with the states outside the trap unpopulated; p to t(0) = v\ ® v 2 <8> p(0), 
where p(0), u 2 = |{0})({0}| 9 i F and v\ = |{0})({0}|^ F are the density matrices of the trapped 
atoms, free molecules and free atoms, respectively. 

Assuming a low loss rate we can expand to second order in k (see, for example, |16|) and 
the equation of motion for the reduced density matrix p = tr F {ptot}, becomes 

-£ = -Jdsjd-xdy /(x,y;t, S )[^ 3 (xt),^( ys )p] + h.c. (2) 

where 

/(x,y;M) = K 2 ({O}|0F(xt)^xt)4(y S )^y S )|{O}). (3) 

Note that the first term in the perturbation, trp {[Hi(s), p to t(t)}} = 0. Following Kagan et 
al.M, we can expand the field operators for the free particles as plane waves modes with the 



energies fi 2 k\/2m and h 2 k\ / 'Am — Ef, for the atom and the molecule, respectively (due to the 
high velocity of the free particles we can neglect the effects of gravity). Since the shortest 
length scale of the trapped atoms is given by the de Broglie wavelength, on the length scale 
of the emitted particles (~ a) , there is no variation of the field of the trapped atoms and we 
can make the replacement ^r(y) ~~ > ^t( x ) in Eq.fl2|). We can then integrate the memory 
function over y which gives rise to the delta function 5(ki + k^), expressing the equal an 
opposite momentum of the molecule and atom. We then perform the integration over k 2 to 
give, 

nt-s) = ^ ? jkldk 1 e-^>-\ (4) 

where F(t — s) = J dyf(x, y; t, s). Note that the x dependence has disappeared. Replacing 
the variable of integration by K = ?>Tik\/A.m in Eq.(f|) and noting that we are only interested 
in frequencies in the range u h — Q to Ub + Q, where E b /Ti 3> Q ^> k, the integral can be 
approximated by 

1_ r +tt AK e -i(K-u b )(t-s) 



F(t - s) w -L / dK e-^-^K*--^ (5) 

where 7 = 3(/t 2 /27r)(4m/3/i) 3 / 2 A /cJb. As £1 is much larger than any system frequencies we 
can write F(t — s) ~ ^Sit — s)/3 and the final Born-Markov master equation for the reduced 
density matrix of the trapped atoms is given by 

d JL = _1_ [H ^ p{t)] + 1 J dx [ 2 ^ 3 (x)p^ t3 (x) - ^ t3 (x)^ 3 (x)p - p^ t3 (x)^ 3 (x)] , (6) 

where I have dropped the subscript on the field operator and returned to the Schrodinger 
picture. This equation is a course-grained equation, valid for time scales of the trapped 
atoms ^> 27r/u; fe and length scales ^> a. By calculating the number decay via Eq.@ we can 
determine that 7 is related to the recombination event rate K3 of Ref. |l ij by K3 = 127. 

This master equation can be used to derive the result of Kagan et al. [§] but it is much 
more general as it can be used to calculate the behavior of any expectation value of the trap 
operators. In this work, however, I am interested in the effect of the three-body decay on 
the quantum state of the condensate. 

For a weakly-interacting gas at zero temperature, the atoms may be assumed to be mostly 
in the condensate. A condensate of n atoms has a wave function that satisfies the time- 



independent Gross-Pitaevskii equation. In the Thomas-Fermi approximation (see [18| and 
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references within) the solution is: gn\ip n (x.)\ 2 ~ /i n — Vt rap (x), for all x where K rap (x) < p n 
and zero elsewhere. The chemical potential p n is determined from the normalization of 
the wavefunction ip n . At higher temperatures the noncondensate particles occupying higher 
energy modes are expected to also act as an environment for the condensate mode. The lack 
of a suitable theory to describe the quantum effects of this interaction (see Ref. |T9|] for an 
attempt) is particularly acute in this case as this interaction is expected to compete with 
the loss from the trap in decohering the state of the condensate. 

Let us assume that the initial condensate states of interest are sharply peaked about the 
mean number N such that \N — n\/N <C 1 and during the short time period of interest 
(the time it takes to lose ~ 3 atoms) the system stays within this linear region. Due to 
the fact that the wavefunction is stable under small changes in number for a > we can 



make the approximations ?/>„(x) « ^jv(x) and \x n ~ /xjv, see Ref.JT^]. The n dependence 
of p n will be discussed latter. Under these approximations we can make the replacement 
ip(x) — > iPn{x-)A in Eq.(H), where A is the annihilation operator for the condensate mode. 
In order to describe longer times where a significant portion of atoms are lost from the 
condensate a more careful approach to defining the condensate mode may be necessary, (see 
Ref. 0). In general, the changing number will give rise to excitations but we can neglect 
these effects for slow decaying numbers over short times. 

With the above assumptions, the master equation for the condensate mode becomes 

^ = 1 hA 3 pti 3 - ti 3 A 3 p - pA^A 3 } , (7) 

where 7 = 7/ dx|?/V(x)| 6 . In Eq.fl7|) I have assumed that the free Hamiltonian has the form 
i?T — TipnA^ A and that we are working in the a frame rotating at the frequency defined by 
pjisf. This equation represents a nonlinear damping of the condensate mode and in the large 
number limit it has an analytical solution. 

Due to the fact that the operator A 3 couples only every third number state, there are 
3x3 independent sub-manifolds of the density matrix corresponding to the number states 
\m;j) = |3m + j), where m is an integer and j = 0, 1 or 2 are the labels for the sub- 
manifolds. The master equation yields an equation for each of the 3x3 the density matrices 
Pjk (j, k — 0, 1 or 2) of the sub- manifolds of the form 

dpjk 



dt ' 1 



2ajPj]~a k — a^cLjpjk — Pjkd\dk (8) 
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where T = 3N 2, y/2. Assuming that the number distribution is sharply peaked at a large 
number the mode operators for each of the three manifolds can be written as 



Y l \/C 3 N + m + IK j) (m+l;j\ + 0(N°), (9) 



lin. 



where C° N = (-2N + 3j — 3)/9, and the sum is over states in the linear region about N. 
The mode operators satisfy [ay,aj,] = 5jji, where I have assumed that only number states 
in the linear region about N are occupied and so it is possible to shift the number states up 
by one. The equations (^) are linear master equations for which their exists a well known 
solution in terms of the coherent states [BJ . 

We can define coherent states for the mode operators aj in the linear region about iV as 



a. 



\<x j ;j)=M'£ , 3 . \m;j), (10) 

im. \](C 3 N + m)\ 

where otj are complex numbers. These states satisfy o,j\ctj] j) « atj\aj;j). They are a 
superposition of every third number state and the number distribution is given by P(n) = 
\{n\atj] j)\ 2 = 5 nt s m+ jQ(m), where, in the large number limit, we can show that 



Q{m) w 



2vr(A^/9) 



9\otj\ 
N 



(11) 



using the relation x\xq x m ^2tcxo e 2x o for x 3> 1 and \x — x \/x 1 p0| . The states 
that we are interested in are peaked at n = N so |<x,-| 2 = iV/9 and from Eq. flTTD we can see 
that the normalization factor in the large number limit is M ~ e _ ' aj ' 2 / 2 . 

The evolution of an arbitrary element of the density matrix in this coherent state repre- 
sentation, \a; j)(j3; k\, is ||: 

(\a;j){i3; k\) t « (a; j) 1 ^ \ae~ Tt - 3 ) (pe~ Tt - k\ , (12) 

where (a;j\/3;j) ~ exp [— ( | ck | 2 + \(3\ 2 — 2a(3*)/2\. For j = k this equation shows that any 
off-diagonal elements in the coherent-state basis, \a;j), rapidly dephase. This means that 
an initial pure state which is a superposition of the states \a;j) will rapidly become a mixed 
state. For example, the superposition of two states with different phases: 

(\a) + |ae**))«a| + (ae**|) — ► |a)(a| + \ae^){ae^\ (13) 
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at the rate ~ |a| 2 I\ This can be seen by noting that {a\ae 1 ^) = exp{— 2i\a\ 2 e 1 ^/ 2 sin(0/2)} 
and therefore the rate at which the off-diagonal terms dephase is ~ |a| 2 r. A number state, 
which is also a superposition of coherent states of different phase, will quickly evolve into a 
mixed state with a binomial number distribution over the number states. For k ^ j, Eq.(|i~2"D 
shows that any state where the coherent amplitude a is not the same in all the sub-manifolds 
will also rapidly dephase. 

The robust states then have the general form \R(a)) = J2j c j\ a ',j), where J2j \ c j\ 2 — 1- 
As the preparation of the initial condensate is via processes that are not likely to preserve 
the independence of the submanifolds, i.e. evaporative cooling, the most relevant states are 
those where Cj = l/v3- in this case, the robust states have the form 

^ a C%+n/3 

\R(a)) = 1 (14) 



where T(x) is the gamma function. These states have the property A 3 \R(a)) = 3Na\R(a)) 
in the linear regime. In general, \R(a)) is not an eigenstate of A but in the large- number 
limit, it is possible to show that (f3\R(a)) ~ e iC " e (P\(3N 2 ay/ 3 ) , where a = \a\e td and \(3) 
is an ordinary coherent state. And so, in the large-number limit, the robust states become 
equivalent to coherent states (except for a phase factor) with the amplitude (3iV 2 a) 1//3 '. 
From Eq . (|i~2j) and |a| 2 ~ N/9, the time scale of the decoherence is 

^ecoh ^ 7^/3, (15) 

where 7 = [(15A) 4 / 5 /14 3 7r 2 ]ii'3(a^aA^) _6//5 . In determining 7 I have used the Thomas-Fermi 
solution of the Gross-Pitaevskii equation for the cylindrically-symmetric harmonic trap: 
Vt(x) = ^mujj_(r 2 + X 2 z), where m is the atomic mass and u± = %/ma\ is the radial 
trap frequency. The decoherence rate is of the order of N times the decay rate for the 
atom number, or, the time it takes to lose ~ 3 atoms from the trap. For example, a 
87 Rb condensate of 3 x 10 5 atoms with K 3 = 2.2 x 10 -28 cm 6 /s and the trap parameters 



uj_ = 2n x 157 Hz and A = 0.075 corresponding to the experiment described in HlO| , yields 
7 w 3.3 x lO^V 1 and r docoh w 4 x 10~ 5 s. 

This decoherence process could be measured directly by observing the rate of decay of 
visibility of the interference of a superposition of states of different phases as a function 
of the lifetime of the superposition. Assuming the condensate is initially in a coherent 



state, a process similar to that described in Ref. [21] could be used to create a superposition: 



coherently transfer half the atoms to another spin state by a | laser pulse; let the phases of 
the two spin states evolve independently and then transfer the atoms back to the original 
spin state with another | pulse. 

To make any conclusions as to the state of the condensate, the above decoherence time 
scale must be compared with that of the collapse of the phase due to condensate-atom 
collisions calculated in Ref.|]T3. Due to atom-atom collisions the chemical potential, fi n , 
is number dependent (this has been neglected in the present work so far) and leads to 
a dephasing of any superposition of number states. Assuming the same trap symmetry 
and the Thomas-Fermi limit of the Gross-Pitaevskii equation as above, the collapse time 
for a number variance of y/N is given by JT3 1 t co \\ ~ (5/u±) (a±/Xa) 2 ^ 5 N 1 ^ 10 . The same 
experimental parameters considered above give a collapse time of r co n w 0.4s, or 4 orders 
of magnitude longer than the decoherence time. In general, r co n/rd C coh oc 0(iV 21 / 10 ) and a 
rough calculation shows that the decoherence process will dominate for > 10 4 atoms. A 
more detailed analysis is necessary but this suggests that collisions will lead a slow phase 
diffusion of a well-defined phase rather than a phase collapse. 

In this work I have presented a realistic mechanism for decoherence at low temperatures 
which supports the assumption that the state of the condensate is a coherent state with a 
well-defined phase. One can also understand the experimental finding that the condensate 
phase is rather robust [2I|; the environment, rather than diffusing the phase helps to keep 
the phase well-defined. However, this work also suggests directions for investigation that go 
beyond the assumption of a well-defined phase. In particular, the observation of macroscopic- 
superpositions of different phase and their decay. In contrast to condensed matter systems 
where there is often a proliferation of possible decoherence mechanisms that cannot be 
isolated from each other, atomic Bose-Einstein condensates are relatively clean systems with 
few possible decoherence mechanisms and are thus a promising test bed for understanding 
decoherence and quantum superpositions in general. Atom loss (single body loss) has been 
shown to have an adverse effect on quantum entanglement between atoms P2|. Therefore 
trap loss (in particular, three-body loss) has important consequences for the possible role of 
Bose-Einstein condensates in quantum information technology. 
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